Generation and Evolution of Spin Entanglement in NRQED 
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57' ABSTRACT 

GO ' 

r**"- . Analysis on the generation of spin entanglement from non-relativistic QED is presented. 

The results of entanglement are obtained with relativistic correction to the leading order of 
(^) 2 . It is shown that to this order the degree of entanglement of a singlet state does not 
change under time evolution whereas the triplet state can change. 
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■ The generation of entanglement states is one of the important problems of quantum 

information science. This is due to the fact that entanglement is the essential source of 
""^ | quantum computing and quantum information processing. An interesting example of en- 

tanglement is the famous Einstein- Podolsky- Rosen (EPR) pairpQ which has been suggested 
to provide secure quantum communication 0. Recently, experimental works have produced 
photons which are in the entangled states and have been used in showing violation of Bell's 
inequality^] and quantum teleportationj^j. However, the generation of entangled states of 
electrons have only been suggested theoretically [S]. By realizing the fact that the long spin 
dephasing time of electrons in semiconductor, the conservation of electron number in the 
non-relativistical limit and the rapid development in the research of quantum dot systems, 
it seems quite promising that electronic entanglement can provide a scalable approach to 
quantum computing. 

Achieving quantum computation requires precise manipulation of the controlled Hamil- 
tonian. In Loss and DiVicenzo's workj^], they proposed to implement quantum gates by 
a time dependent Heisenberg spin Hamiltonian. More recently, generation of EPR states 
from two-electron mixed state has also been proposed in quantum dots system j^j. As all 
effective interactions in condensed matter system are basically electromagnetic, it is there- 
fore interesting to address the entanglement problem within the framework of Quantum 
Electrodynamics(QED). This approach not only has the advantage of relativistic invariance 
but also takes care of the problem of indistinguishable property of electrons. In a recent 
paper, Pachos and Solano [7j have discussed the problem of the relativistic invariance of en- 
tanglement by using QED. They have also claimed to obtain the Heisenberg Hamiltonian by 
considering the lowest-order scattering process of two electrons. Since magnetic interactions 
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are of relativistic nature, it is important to consider all the relativistic effects to the same 
order. Due to the fact that electrons are fermions, the antisymmetric nature of the state 
vector leads to interesting results which do not arise for non-identical particles. Moreover, 
there is also another aspect of entanglement which needs to be addressed. Namely, the 
effects of interaction on the evolution of entanglement states should be investigated. The 
problem of entanglement and identical particles has also been addressed by Omar et al[H], 
however they concentrate on transferring entanglement from the internal to the spatial de- 
grees of freedom without considering interaction. In this work, we analyze the generation 
and evolution of entangled electron pair within the non-relativistic(NR) expansion by in- 
cluding the interaction effects to the order of (^) 2 and using properly antisymmetric wave 
function. Consistency then requires the QED corrections to 0(a 2 ) with a being the fine 
structure constant. Therefore, for completeness, the two-photons exchanged processes are 
also discussed in this work. 

Due to the non-relativistic approach of this work, the state of the two-electron system 
can be expressed by a wave function. Furthermore, with electrons being identical fermions, 
the total wave function of the electron system is required to be totally antisymmetric. In 
discussing scattering process the total wave function of two electrons VP can be expressed 
by the direct product of spatial wave function ijj and spin wave function \ '■ ^ = V* ® X- 
For the moment it is appropriate to introduce the concept of Schmidts decomposition on 
the spin wave function such that the entanglement analysis can be performed effectively. 

The wave function of any bipartite system is expressed as a double sum. The Schmidts 
decomposition procedure asserts that the double sum can be expressed as a single sum by 
local unitary transformation. For example, the spin state vector of two spin-1/2 particles 
is: 

\x) = ^2C nm \( n )A\Urn}B, (1) 
nm 

where {|Ci)> |C2)}a and {|<^i}, \^2)}b are two sets of orthonormal basis which belong to the 
spin Hilbert space of the bipartite system A and B. By applying the Schmidts decomposition 
procedure |x) can be reduced as : 

Ix) = a,i\rn) A \£i) B + a 2 \m)A\i2)B- (2) 

Here, |??2}}i4 and {|£i), |£2)}b are the two possible incomplete orthonormal basis of 

the Hilbert space, where {01,02} are Schmidts coefficients which satisfy the normalization 
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condition of = 1. This decomposition procedure is particularly useful for expressing 
i=i 

entangled state. For this two-spin system, if all aj are non-vanishing then the system is 
spin entangled which means the state vector \x) can not be expressed as a direct product 
W)a\C')b- By picking a special direction as z-axis, one can define a conventional basis {|| 
), ||)} along this axis such that any spin state vector can be expanded as linear combination 
of these basis vectors. Denoted a symmetrized and entangled spin vector as \xs)> the general 
form written in the conventional basis is: 

lxf> = C|TT>+£>(IU>+liT>) + G|U> with^VCG. (3) 
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By tuning the coefficients C, D and G, three of the Bell states, = -i=(| jj)± III)) an d 
|^ + ) = -^(|T|)+ ||T)) can be obtained from Eq.(3). If D 2 = CE then the state is not 
entangled, that is to say by local unitary transformation the state can be written as a direct 
product. On the other hand, for the antisymmetric and entangled spin vector \Xa) one can 
easily prove that the state is the remaining Bell state l^ - ), 

\Xa) = |*~) = ^(lti>- lit)), (4) 

It is also noted that if the spin degree of the system is not entangled, then one of the 
Schmidts coefficients must vanish, and the other is equal to one. 

By properly anti-symmetrized the total wave function and choosing the Coulomb gauge, 
the interaction potentialjlUj to the v 2 /c 2 is: 

U = U c + Usl + U ss , (5) 

where Uc, Usl and Ugs, denotes the Coulomb interaction with relativistic correction, 
spin-orbital interaction, and spin-spin interaction respectively. The form of these three 
interactions are as following: 

"o = - - g^(r) - ^{ Pl ■ p, + r • (r • 2 Pl)P2 } (6) 

r 4crm z 2m z c z r r z 

e 2 h 

U ls = - 4m 2 c 2 r3 { r x (Pi " Pa) • (*i + ^2)} (7) 
e 2 h 2 rr 87r . 1.. , 3(<Ti • r)(cr 2 • r). 

4m z r 3 r J r° 

It is noted that all the correction terms are of the same order which is v 2 /c 2 correction to the 
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Coulomb potential By applying these interacting potentials in the Born approximation, 
a systematic analysis on entanglement can be proceeded. However, it is not consistence if one 
stops at the first Born approximation of order a. As well known in atomic physics problem, 
the energy correction due to the relativistic effect is the same order as a 2 . Furthermore, 
for electrons in metal, the Fermi velocity vf is of the order of 10 6 m/s and hence the ratio 
vf/c is the same order as a. Therefore, consistency requires the consideration of the second 
Born term which is 0(a 2 ). This higher order analysis will be discussed later. 

In order to address the generation and evolution of entangled states of two electrons, 
the effects of the spatial wave function must be considered. According to the principle of 
quantum mechanics on identical particles, the total wave function is either symmetric or 
antisymmetric depending on the nature of the system. Therefore any sensible discussions 
on entanglement must take into account the spatial property. This is in contrast to the 
conventional discussion on entanglement where only the spin states are involved. With this 
in mind the problem can be posed in the following way. If the system is prepared without 
spin entanglement, then after scattered by the above interaction potentials, is it possible 
to generate a spin entangled state? The more interesting question is what type of spin 
entanglement state is robust during the scattering process? These questions are addressed 
in the next section. Let us consider the two electron scattering process. The initial wave 
function can be prepared as the simultaneous eigenfunction of the Hamiltonian and total 
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momentum operators. The three possible different forms of the wave function are: 



l* (1) > = hMPi,P2))®|x?> 

l* (2) > = l^(Pi,P 2 ))®|xf> (9) 
l^ (3) ) = l^(pi,P2))^|xi)- 

Here, \ips (Pi;P 2 )) and |Vvl(Pi> P 2 )) represent respectively the symmetric and antisymmet- 
ric wave functions of the system with momenta pi and p2- (Such initial states in principle 
can be created by doubly ionizing the two electrons in the ground state of the helium atom.) 
The explicit form of these functions \ips) and are : 

l^(P!,P2)) = — i^{ e i(Pi^+P- x =) +ei(P- x ^^)} 

lVu(pi,P2)) = (27r ^3/ 2 i ei(Pl ' Xl+P2 X2) " e^P-^+P 1 ^)}. 

The form of |xf ) and \xa) are gi yen by Eqs. (3) and (4) respectively. The state \xs) 
representing the no entangled spin state is : 

|xf) = C|TT)+f(IU>+liT» + G|U) with d 2 = cg. (10) 

Let the initial and final states of the system be denoted respectively as l^f ) and \^f^) 
where k represents any one of the states of Eq.(9). Due to the fact that some of the terms of 
the interaction potential U, such as Uc, Uls and {a± • 02}; are commuting with total spin 
S = Si +S2, therefore the scattering process through these terms is spin angular momentum 
conserving. Even though the last term of Eq.(8) does not commute with the total spin S, 
however when acting on either symmetric or antisymmetric spin state, the resulting state 
retains the same symmetry property. This is due to the fact that this term is symmetric. 
Explicitly the form of the resulting states are: 

(<7 1 T)( < 7 a .r)|xA) = -^|XA) (11) 



(ax • r)(a 2 • r)| X s) = wi |TT) + w 2 (|U) + lit)) + ^3 III) (12) 

where U{ are some spatial factors. Then it is obvious that the scattering between symmetric 
and antisymmetric spin state by this term produces zero transition amplitude. As a result 
the total scattering amplitude of the initial state |^^) is : 

This result shows that as the initial state being no spin entanglement, the probability 

(3) 

for creating a spin singlet entangled state ) is zero. This result also holds for the 

(2) 

initial state being ), which is an entangled spin triplet state. On the other hand, 

(3) 

by taking the spin singlet state ) as initial state, one can make use of the complex 
conjugate of the above results to conclude that the final state can only be the same spin 
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Figure 1. (a)The ladder and (b)the crossed diagrams of two-photon exchange processes. 



entangled state. Furthermore, the results also imply if we prepare an initial spin state as 
either \xs) or \Xs)> the interaction potential U can change the degree of spin entanglement. 
For example, the spin triplet entangled state \xs) can scatter into a no spin entangled 
state \Xs)- Therefore one has an important conclusion that, by taking QED into account, 
the degree of entanglement is not an invariant concept during time evolution. From above 
discussions, one may conclude that, due to the stability of entanglement of the spin-singlet 
state, implementation of spin-singlet state in quantum computation is more practical. 

In order to justify the above conclusion, it is necessary to establish the fact that the 
high-order corrections do not spoil the result. Furthermore, there is also another important 
issue which requires at least to include the next order correction. This was mentioned earlier 
that all correction terms are of the same order 0(v 2 /c 2 ) correction to the Coulomb potential 

Since the magnitude of v 2 /c 2 is the same as a 2 , the second Born approximation by the 
interaction potential U should be included in the discussion. 

A complete a 2 correction to the scattering amplitude must include also the radiative 
correction to the coupling constant a and the electron mass. However, since this work is only 
considering the non-relativistic expansion, these radiative effects do not alter the results of 
scattering and will be neglected. The second Born approximation can be obtained directly 
from calculating the corresponding Feynman diagrams which are the ladder and crossed 
diagrams given in figure 1. This approach is quite elaborated for the result of this section. 
Fortunately, there exists an effective way of calculating the results within the validity of the 
non-relativistic expansion. It is known that the a 2 -correction can be obtained by a second 
order potential : 

UW = VP + VP, (14) 

(2) (2) 

where V£ and V£ are calculated from the ladder and crossed diagrams respectively. Even 
though the form of these potentials are gauge dependent, it has been shown that the total 
scattering amplitudes are gauge independent |llj. Therefore one can chooses a convenient 
gauge to calculate the amplitude. In the present study the Coulomb gauge is more suitable, 

(2) 

as will be shown latter, the leading relativistic correction of the vanishes in the Coulomb 
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gauge. 

To ease the discussion, the notation for the wave function is modified slightly. All the 
initial states of two electrons are written as \ab) and the final states as \cd). The lowest 
order scattering amplitude (cd\S^\ab) can be identified with a potential scattering with 
the interaction potential U given by Eq.(5) : 

(cd\S m \ab) = -2Tri5(e a +e b -e c -e d )(cd\U 12 \ab), (15) 

where e a denotes the energy eigenvalue of the free electron of state a and U12 is the potential 
of vertices 1 and 2 . The two-photon exchange contributions to the order of interest, namely 
0(a 2 ), can be expressed as the effective potentials which are given by : 

(cd\ V P\ab) = -—y[ dz {C3d ; lU ^ £a ~ £c ~ ^^(r^liW^N^) 
2-ki J [£ a — e r — z + ir]E r ) (£fe — e s + z + vqe s ) 



and similarly for the crossed diagram 

(cd\V {2) \ab) = 1 V f dz {c-iS 4 \U M {£ a - e c - z^r^jnsd^P^jz^a^) _ 

X 27ri TfJ (e a -e r -z + irie r )(e d -e s + z + irie 3 ) 

Here, Uij(cq) is the Fourier transform of the potential with q being the three dimensional 
momentum transfer. As shown in Lindgren's work^J, the result of Eq.(16) can be simplified 

as : 

/,, T/ (2)| ,v _ 1 {{cd\r^\rs){rs\V BM \ab) + {cd\V BM \rs){rs\r^\ab)} 

{cd\V L H-^L {£a + £h _ £r _ £g) ( 18 ) 



Here, V BM is the generalized potential of Brown and Mittlemanjl2j 

{rs\V BM \ab) = \{rs\U{e a - e r ) + U(e b - e s )\ab). (19) 

Since what one needs from this calculation is the order a 2 correction, it is obvious that 
one should only keep the Coulomb potential e 2 /r of U in the above expression for V BM . 

However, due to fact that the Coulomb potential does not flip the spin of electrons, hence 

(2) 

V£ does not scatter the spin-singlet state into other spin configuration. It is now important 

(2) 

to show that also preserves the spin-singlet state. 

For the crossed diagram, one can also simplify Eq.(17) and the result is: 

(d\V {2) \ b)= 1 V {(cs\r^\rb)(rd\V di ff\as) + (cs\V di ff\rb)(rd\r£\as)} 

{ ' X ' ; " 2^ ( £a+£s - £d - £r ) ^ > 

, where V d% ^ is the difference potential 

(rs\V di ^\ab) = -(rs\U(e a - e r ) - U(e b - e s )\ab). (21) 

Thus, within the same approximation, one should retain only the Coulomb potential 
and therefore a vanishing V dl ^ is obtained. As a result, to the leading relativistic order, 
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does not contribute to the scattering process. This completes the proof of the stability 
of the spin-singlet entanglement. 

Recently, people have tried to control electron spin as qubit directly in some condensed 
matter system, such as electrons in quantum dots. It is also claimed the interaction of such 
system can be described by the Heisenberg Hamiltonian H = J(jb)o\ a?,. It is easy to check 
that the Bell states are stable during time evolution with this Hamiltonian. However, in 
this approach, spins are the only variables of the system which, in general, can not properly 
describe electrons. As pointed out earlier, it is inadequate by neglecting the spatial wave 
function which is involved in the formalism of identical particles. The effects of identical 
particle on entanglement can be seen from previous discussion that the entangled spin-triplet 
states, such as Bell states, can evolve into no spin entangled state. As a result the concept of 
entanglement for spin-triplet state is not practical. This is in contrast to the case described 
solely by the Heisenberg Hamiltonian. On the other hand, as shown also in this work, 
the spin-singlet state which is entangled by nature is stable to a high degree of accuracy. 
Therefore, for any practical implementation of entangled state with spin, the spin-singlet 
state should deserve more attention. Furthermore, the result obtained here by considering 
the antisymmetric wave function, implies that the spin-singlet entangled state can not be 
generated from the other spin configuration regardless of its state of entanglement . 
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